The microscopic spin-phonon coupling constants in CuGeC>3 



Ralph Werner and Claudius Gros 
Institut fur Physik, Universitat Dortmund, D-44^21 Dortmund, Germany. 

Markus Braden 

Forschungszentrum Karlsruhe, Institut fur Nukleare Festkorperphysik, Postfach 3640, D-76021 Karlsruhe, Germany 
and Laboratoire Leon Brillouin (CEA-CNRS), Centre d'Etudes Nucleaires de Saclay, F-91191 Gif-sur-Yvette Cedex, France. 

(February 1, 2008) 



On ■ 
On ; 

or 

<■ 

(N ; 

> 
oo 

CO 

o 
o 

00 



c 

o 
o 



Preprint. Typeset using REVTj^X 

Using random phase approximation results, mean-field 
theory, and refined data for the polarization vectors we deter- 
mine the coupling constants of the four Peierls-active phonon 
modes to the spin chains. We then derive the values of the 
coupling of the spin system to the linear ionic displacements, 
the bond lengths and the angles between bonds. Our values 
are consistent with microscopic theories and various exper- 
imental results. We discuss the applicability of static ap- 
proaches to the spin-phonon coupling. The c-axis anomaly 
of the thermal expansion is explained. We give the values of 
the coupling constants in an effective one-dimensional Hamil- 
tonian. 



The microscopic three-dimensional (3D) Hamiltonian 
we have to consider consists of three parts: 



I. INTRODUCTION 

With the discovery of the first inorganic spin-Peierls 
compound CuGe03,El it has become possible to investi- 
gate the physics of the spin-Peierls transition in quasi- 
one-dimensional spin chains with high precision. As a 
consequence, it has attracted much attention both in ex- 
periment and theogu. Yet, for a long time the absence of a 
soft phonon modeou has been puzzling, since the behav- 
ior was believed not to be consistent with the standard 
approach to spin-Peierls transitions by Cross and Fisher .□ 
The frequencies of the Peierls-active phonon modes be- 
ing of the order of the magnetic exchange^ it has been 
argued that the Cross and Fisher approach is not-appli- 
cable because of the nonadiabagity of the phonons.Q Only 
recently we were able to showu that the random-phase- 
approximation (RPA) approach by Cross and Fisher in- 
deed is consistent with the hardening of the Peierls-active 
phonon modes .El 

It is now tempting to combine the RPA results with thft 
detailed data on the phonons acquired by Braden et alia 
Treating the lattice with the standard harmonic theory 
and including the spin-phonon coupling mean-field like, 
we calculate the microscopic coupling constants between 
the lattice and the spin chains. It is then possible to pre- 
dict the effect of structural changes on the spin system, 
namely the antiferromagnetic exchange— I. The latter has 
been suiuect to various experimentalfJlilJ and theoretical 
studies-EliZl 



H = H X + H„ 



H sp . 



The Heisenberg spin Hamiltonian 

H s = J ^ Si • Si + £ + J 2 ^ Si • Si + 2£ 



(1) 



(2) 



with the exchange integrals J and J2 between nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) Cu d 
orbitals, respectively. 

Further we distinguish the phonon part 

ff P = E%^+ E K'£,„yrZ tV ri y , (3) 



describing the lattice vibrations in the harmonic approx- 
imation, where r n „ = (rj „, v , „) are the deviations 
from the ionic equilibrium positions. 

Finally the spin-phonon coupling term reads 



H, 



sp 



1 



A Ji,i+z Si ■ Si + £ + A J\j+2z Si • S 



1 + 22- 



(4) 

The energy scale A Ji,i+£ is a function of the variation of 
the magnetic exchange integral with the atomic displace- 
ments gf v — dJ/dr" v to be discussed in Sec. III. The 
NNN term AJi 1+22 is a function of dJ-xldr^ v . 

The indices used are n = (n x ,n y ,n z ) G Z 3 running 
over all unit cells of the three-dimensional crystal, the 
Cu-site index 1 = (l x , L, l z ) £ Z 3 (two Cu sites per unit 
cell), and the unit vectors x = (1,0,0), y — (0,1,0), 
and z = (0, 0, 1) to nearest-neighbor unit cells in the 
corresponding direction. The index v labels the 10 atoms 
within a unit cell as shown in Fig. j] and a € {x, y, z} is 
the vectorial component of the indexed quantity in the 
respective three-dimensional space. 

In Sec. |n] we briefly summarize the diagonalization of 
the phonon Hamiltonian (^) followed by the discussion of 
the symmetry of the four Peierls-active phonon modes, 
including refined data for their polarization vectors. Us- 
ing these symmetries we transform in Sec. DTD the micro- 
scopic spin-phonon coupling Hamiltonian (^) to normal 
coordinates in reciprocal space. This procedure yields 
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relations between the different linear, angular, and nor- 
mal mode coupling constants. Using RPA results and 
mean-field theory in Sec. [v| we obtain numerical values 
for the normal-mode coupling constants which then can 
be converted to the real-space coupling constants. The 
resulting dependence of the magnetic exch ange on static 
distortion of the lattice is discussed in Sec. VII and com- 
pared with values from the literature. Finally we derive 
an effective one-dimensional model to give coupling con- 
stants consistent with frequently applied theoretical ap- 
proaches. The consistency of the different results gives a 
a posteriori justification of the mean-field approach. 



II. PEIERLS-ACTIVE PHONON MODES 




In the standard treatment of harmonic lattice dynam- 
ics, the initial problem of 3 • N ■ Ni on degrees of free- 
dom (A number of unit cells, Aj on number of ions in 
the unit cell) is transformed into reciprocal space by a 
Fourier transformation, where-JV wave vectors fulfill the 
periodic boundary condition.Ej For any fixed wave vec- 
tor one obtains a 3 • Ai on -dimensional problem which may 
be diagonalized, resulting in a set of 3 • A; on eigenmodes 
labeled by A S {I, . . . , 3Ai on }. For that purpose, the 
displacement and momentum operators are decomposed 
into eigenmode contributions introducing normal coordi- 
nates Q and conjugated momenta P: 

v q a 

^ q a 

The vectors R n designate the coordinates of the unit- 
cell origins, rriy is the mass of the z4;h atom and e v are 
polarization vectors. Note that we use a nonstandard 
definition for R n and e„(A, q) which will simplify the in- 
terpretation of the polarization vectors at high symmetry 
points in the Brillouin zone. Further transformation to 
boson creation and annihilation operator representation 
via 

QA,q=^(^ q + 6 A , q ), (7) 

^ = V^( 6 U-<)' (8) 

yields the Hamiltonian usually used in the theoretical 
treatment of the lattice vibrations 

q,A V J 

With the experimentally determined phonon modes 
f^A.q and shell-model calculations, it is possible to de- 
termine the components of the polarization vectors 



FIG. 1. Projection of the unit cell of CuGeOa on the x-y 
plane. The oxygen atoms are distinguished into 01, 02a, and 
02b, the atoms of the second formula unit are labeled with 
a prime. Each unit cell contains two Cu chains in z direction 
(positive 2 direction is into the plane). The broken lines show 
the reduced unit cell introduced in Sec. n is the index for 
the whole cells, 1 indexes the reduced cells. 

e„(A,q).i At the wave vector of the Peierls instability 
qo = (7r/a, 0, ir/c), four of the 30 modes correspond to 
the irreducible representation with the symmetry of the 
lattice distortion in the spin-Peierls phase, Tj in the no- 
tation of Rcf. p|. a = 4.8 A, b = 8.5 A, and c = 2.9 A 
are the lengths of the unit cell in x-, y-, and z-direction, 
respectively. 

Adapting the lattice-dynamical model presented in 
Ref. ||, a special effort was made for the description of 
the spin-Peierls relevant modes by the introduction of ad- 
ditional force constants; details of the lattice-dynamical 
study will be given elsewhereEj The T^" modes are char- 
acterized by displacements of the Cu ions along c, of the 
02 sites along a and 6, and of the Ge ions along b. The 
polarization pattern of the four modes as obtained by the 
shell model are represented in Fig. ^ and given together 
with their frequencies in Table Q. The highest T^ mode 
corresponds to a Ge-0 bond stretching vibration thereby 
explaining its elevated frequency. The three T^~ modes at 
lower energies posses a common element which consists 
in the rotation of the 02-02 edges of the Cu04 plaque- 
ttes in the x-y plane around the c axis. However, only for 
the lowest mode does this twisting of the Cu02 ribbons 
describe the main character of the polarization pattern. 
The modes at II and at 6.5 THz show, in addition, a 
modulation of the lengths of the 02-02 edges and a Cu 
shift parallel c. For the 11 THz mode the displacements 
of the Cu ions modulate the 02-Cu bond distance. The 
6.5 THz mode is characterized by a strong modulation of 
the Cu-02-Cu bond angle (see Fig. ||) which is essential 
for the magnetic interaction. 

For later use we define the matrix M with the elements 
given by e"(A, qo)/ \jvn v extracted from Table Q: 
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24.6Thz 



ll.lThz 



6.53THz 



3.12THz 



FIG. 2. Geometry of the Tf eigenmodes as given by the 
polarization patterns in Table \A The shaded areas are the 
CUO4 plaquettes which form the Cu chains in z direction. 
The Cu atoms are in the center of each plaquette, the corners 
are formed by 02 ions. The Ol atoms are represented by the 
open circles with the Ge ions in between them. (Compare 
with the x-y projection given in Fig. jjj.) Note that the 02 
elongations are in the x-y plane, the (small) Cu displacements 
are along the z Axis while the Ge displacements are along y. 
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0.0095 
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-0.4790 


0.7412 


-0.0888 




-0.4330 


-0.5325 


-0.3698 


-0.2605 


weg 2b (A)/V^o2 


-0.6212 


0.6581 


0.3382 


-0.7932 




-0.8620 


0.1339 


0.2021 


0.8723 



TABLE I. Frequencies and polarization of the Peierls-ac- 
tive T^" phonon modes at room temperature. The global pref- 



actor is given by u 
is e5(A,qo) = e?(A). 



:.26 ± 0.02) x 10" ib kg. The notation 



M = 



/ 0.03 -1.67 2.58 -0.31 \ 

-1.51 -1.85 -1.29 -0.91 

-2.16 2.29 1.18 -2.76 

V -3.00 0.47 0.70 3.04 ) 
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(10) 



The static distortion in thc_apin-Pcicrls phase at T = 4 K 
also has been determined £j We define a corresponding 
four-dimensional vector: 



(r) 



T=4K 



( (rL) \ 

( r 02b) 



10" 



/ 0.57 \ 
0.08 

-0.95 
V -0.65 J 



A. (ii) 



III. SPIN-PHONON COUPLING TERM 

In the spin-phonon coupling term (|J) we focus on the 
NN part for reasons that will become obvious at the end 
of the section. We include the relevant displacements of 
the ions directly involved in the Cu-Cu supcrexchangc 
path determining J and only those coupling constants 
where the ions actually show displacements in the Peierls- 
active modes. The apex Ol atoms are not displaced by 
those modes at the appropriate wave vector qo. We have 
to consider two copper ions in adjacent unit cells along 
the c direction, two germanium sites, and two oxygen 
atoms surrounding a Cu-Cu bond. The notation intro- 
duced is shown in Fig. [l], the two 02 oxygen atoms per 
formula unit are denoted 02a and 02b. There are two 
formula units per unit cell which we distinguish by a 
prime. 

The relevant coupling constants for the linear atomic 
elongations are shown in Table [n| The effective spin- 
phonon coupling Hamiltonian is 



H 



NN . 



E 



Cu,n+2 



8 )-pL(tL', 



3 



- — r 
-y 



y 

02b,n 



+ E 



- a v (r y 

5c u ( r Cu',n — r Cu',n+l) — 5gb ( r G 

"h 5o2a ( f *02a,n f *02a',n) 
~ i/02a I ' 02a,n ' Q2a' , n J 



Gc,n Gc',n 



• S' a . (12) 



The two sums correspond to the two Cu chains running 
through each unit cell. 

The symmetry of the Hamiltonian (|l^) allows for some 
simplifications. First of all we use the equivalence of 
coupling to the 02a and 02b displacements (see Fig. [I]): 



5 2 
f 2 



5o2b 5o2a- 



(13) 
(14) 

From the symmetry of the Tj modes (see Fig. |J) we see 
that the 02-y components are in phase, i.e., 7 , 02an = 
^o2b,n = r o2,nj while the x components exhibit an an- 
tiphase shift: -r£ 2a n = --r% 2b n = rg 2>n . 

As indicated in Fig. fi] we then cut the unit cell along 
the y axis in half separating the ions labeled "prime" 
from those without a label. 



S' • S 



n+5 



r a 

i+y/2 ' S n+ y/ 2 +z- 



(15) 
(16) 



The change of sign of the coordinates accounts for the 
antiphase elongation of the two types of ions in the 
Peierls-active modes at the wave vector of the instability 
q = qo = (7r/a, 0, n/c). Resummation n — > 1 over all the 
new cells, i.e., twice as many with a new cell length 6/2 
in the y direction, yields 



H. 



NN . 
sp 



=E 



1 L 







+ 3S2 


J"o2, 1— y 


1- r 03 ,l) 








+ 9 V 02 


(r v 

\ 02,1-3/ 




Si +S . (17) 






dJ 






dJ' 




9r Cu,„ 




° Cu,n+z 




gL 




dJ 




dJ 


dJ' 




*■&.,» 




Ge',n-$ 


Go,n 


x.y 




dJ' 




SJ' 


aj' 




dr x - y 
02a, 




02a',n 


dr x ' y , 

02a' ,n-y 


x,y 
9o2h 




dJ 




SJ 


8J 




Q2b, 




Q2b',n 


dr x,y , 

02b',n-y 



TABLE II. Definition of the coupling constants for linear 
atomic elongations. The two Cu chains running through each 
unit cell are distinguished by a prime (see Fig. Q). J' is the 
magnetic coupling constant along the Cu' chains. 



The overall change of sign in the first sum with respect 
to the second in Eq. (|12| ) has been incorporated in the 
phase factor e mly . This change of sign translates into 
the antiphase shift of the spin-Peierls ordering between 
neighboring Cu chains in the y-direction. 

Now we substitute the displacements r" t with the q- 
space normal coordinates (^|). For clarity we introduce 
the abbreviation 



(18) 



for the Fourier transform of the nearest-neighbor spin- 
spin correlation operator: 



H 



NN 
sp 



^E s -« E 



2fi 



A.q 



3A(q)QA, q - (19) 



q a 

Here the effective normal mode coupling constant 

e Cu (A,q) 



2n 



A,q 



g x (q) := (1 - e*« c ) 



+ (e-W + l) [gv. 



eL(A,q) 



a: eg 2 (A,q) , ny e^ 2 (A,q) 



.9" 



(20) 



1 JUi > 

was introduced@ 

The next step is to transform the normal coordinates to 
boson creation and annihilation operator representation 
via Eq. (@): 

E ™ = ^f=L ^ E *A(q) (*!,-, + *>,,)■ (21) 



This is the representation usually used in theoretical ap- 
proaches to spin-phonon coupling.!] Since the polariza- 
tion vectors are known for qo (Table ||) , Eq. (^0|) defines 
the relation between the coupling to the linear atomic 
deviations g" (Table [il]) and the normal mode coupling 
constants g\ (q) . 

The NNN exchange term J2 Si • Si+2,2 leads to a mag- 
netoelastic coupling equivalent to the one for the NN ex- 
change shown in Eq. (|l_7j) . Including all ionic linear elon- 
gations contributing to the Cu-02-02-Cu NNN superex- 
change path the prefactors in the resulting reciprocal- 
space coupling constants — compare Eq. ( pfj| ) — then 
are (1 - e 2iq * c ) for the Cu part and (1 + e l « c ) for the 
other ions. The coupling of the J 2 term vanishes at the 
wave vector of the instability q = qo = (ft/a, 0, tt/c) and 
does thus not directly influence the transition. Within 
the present work we can disregard this contribution. 



IV. COUPLING TO BOND ANGLES AND 
LENGTHS 

The two lower Peierls-active modes essentially vary the 
angles r] R = Z[Cu-02k-Cu] and S K = Z[02«/-02K-Ge]. 
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V 
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\d6 K /dr z Cu \ 


0.11 


71 

A 







0.76 





\d6 K /dr y Ga \ 







0.11 


TV 

A 





0.82 


\d9 K /drl 2 \ 


0.21 


71 

A 


0.29 


7T 

A 


0.54 


0.57 


|W"/Srg a | 


0.14 


7T 

A 


0.32 


7r 
A 


0.36 


0.82 


Oo 


0.55 


7T 


0.89 


7T 


1.93 A 


1.73 A 



TABLE III. Linear coefficients of the expansion of the an- 
gles and bond lengths as a function of the linear atomic elon- 
gations, as defined in Eq. (^). The variables k G {a, a , b, b'} 
and 8 K £ {77", 5 K , d£ u , d GB } are introduced in the text. The 
last line holds the experimental equilibrium angles and bond 
lengths (Ref. [H]) [9 K =9 + A8 K , A8 K is defined in Eq. (§|)]. 



Together with the bond lengths 



dL = 



Cu — 02k and 



Ge — 02k they represent the natural set of coor- 
dinates of the lattice vibrations in the irreducible group 
of the Tj modes. The index k £ {a,a',b,b'} was intro- 
duced to label the position on the different oxygen atoms 
in the unit cell. 

Introducing the variable 9 K € {rj K , 5 K , d£ u) dJ5 e } we can 
write 

^AO- = X^n n X^- ^e " 
^6z, 



A J, 



1,1+2 



{9} n {8} 



dr a 



(22) 



Here we defined the coupling constants ge = (dJ)/(86*), 
which are independent of k. For reasons of translational 
invariance we can drop the unit-cell index n. The lin- 
ear coefficients of the Taylor expansions {d0 K )/(r") at 
different positions k in the unit cell all yield the same 
numerical coefficients but with varying signs. The abso- 



lute values of the coefficients are given in Table III. 

Considering all the relevant bonds and angles and us- 
ing the decomposition (|22) we can set up a spin-phonon 
Hamiltonian similar to Eq. (|l^) in the previous section. 
By a simple comparison of the coefficients we obtain the 
transformation matrix between the angular and bond- 
length coupling constants and the linear atomic deviation 
coupling constants. 



9c 

i/Go 



\ / -0-22 I 



9o2 

\glj 









0.21 f 

A 

V -0-14 f 



96 



-1.52 
-0.11 f 0.82 
0.29 f -1.08 0.57 
0.32 f 0.72 -0.82 / \ ffl / 

(23) 



9t 



Together with Eq. (|2(j) we now can determine all coupling 
constants if any four of the them are known. 



V. NORMAL MODE COUPLING CONSTANTS 

We now numerically determine the four normal mode 
coupling constants. .-We have shown the RPA approach 
by Cross and FisherB for the Hamiltonian 



H — J Si ■ Si + £ + J 2 2_j ^ 1 ' ^1+21 
1 1 

q,A V 

^E S -« E 5A(q) + 6 



A.q 



(24) 



to satisfactorily describe the dynamics of the Peierls- 
active Phonon modes.ll It consists of the Heisenberg chain 
(^), the harmonic phonon part (^) and the spin-phonon 
coupling term (|2l|) all discussed above. We have given 
an expression fot the critical temperature of the spin- 
Peierls transition!!! which we generalize to the four Peierls- 
active modes_we have to consider herein qo = 
g x (q )^g x )H 



25! 



2sf 



2 5 1 



fc B Tsp= T^ + T^ + T^ + T^ XO (25) 



Ml-, 



The factor xo ~ 0.5 is a contribution of the static spin- 
polarization function at the appropriate wave vector. Its 
value is controversial and we have adopted a mean of the 
proposed values. Please refer to the discussion in Sec. IX 
for the details. 



A. Ginzburg criterion 

The Ginzburg criterion gives an estimate of the tem- 
perature range of the critical region in which fluctuations 
suppress the applicability of mean-field approaches (or 
RPA). It is obtained through comparing the theoretical 
correction of Gaussian fluctuations to the specific heat 



P,0 



T 2 

1 SP 



abc 



16tt (T - Tgp) 2 Ubtc 



(26) 



with the experimental jump in the specific heat at the 
transition.^ 2 ] The correlation lengths £„, and £ c along 
the respective crystallographic axes can be obtained from 
fits to the diffuse x-ray data from Schoeffel et al£3 

£ a «0.50a [(T-T SP )/Tsp]-*, 
6 « 0.65 b [(T-T SP )/T SP ]-i , 
£ c « 3.06 c [(T-T SP )/T SP ]-i . 



■■heat jump has been determined by Lasjau- 
i to be AC, 



The specif 
nias et all 
cell volume. Requiring C p 
Ginzburg criterion to be 



cxp 



= 0.73fce at Tgp per unit- 
Cp t o <C AC cxp we find the 



(T-Tsp) > 0.03 T SP = 0.4 K. 



(27) 



In accordance with the meaia-field approach to the sus- 
ceptibility by Kliimper et al!EB we conclude that beyond 
a region of 3-4 K around Tsp the mean-field theory is 
reliable. 
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B. Mean-field approach 

The transition temperature given for CuGe03 with 
14.1 K, one parameter is fixed through Eq. (p5|). As we 
shall discuss now, the others can be estimated from the 
polarization vectors of the Peierls-active phonon modes 
and the static distortion lb-, the dimerized phase at 4 K 
also given by Braden et aZli3 For the fixed wave vector of 
the Peierls instability q we can derive from expressions 
(|J) and (0) a relation between a static lattice distortion 
(r") and the expectation values of the displacement of 
the eigenmodes out of the harmonic equilibrium (i>A,q ) : 



(rqo,,) = e°(A,qp) 



2h 



N 



0> 



v - ■ ■ y - £ A,qo 

Introducing the canonical transformation 

1 0x(q) a 



(Kqo)- (28) 



?A,q 



3A,q 



A,q 



(29) 



for the Bose annihilation and creation operators b\ „, 
where S q was defined by Eq. (|l8|), the Hamiltonian j24] ) 
decouples into 



\g>M: 

Aq 



E^q^A.q ' , 
q,A X 



(30) 



The operators &A,q do not satisfy Bose commutation re- 
lations and since [S-q, b X q \- ^ the solution of this 
Hamiltonian is not at all evident. But in a mean-field- 
like approach we can assume (b x ) = so that from Eq. 
(p9|) follows 



1 5a (q) 



N M} A ,q 



(31) 



The mean-field ansatz is reasonable here, since we are in- 
terested in temperatures of 4 K which is far from the crit-. 
ical region and the dimerization is as good as saturated.^ 



C. Values 



With Eqs. (E8|) and (|31j) we are left with a set of linear 
equations. The values of the frequencies £l\ are given in 
Table H, the polarization vectors e"(A, qo) enter via the 
matrix M defined in Eq. ([To|), and {t)t=4k is given in 
Eq. P). 



(r) 



T=4K 




(32) 



gi/ks 


g2/k B 


93 Ab 


34 Ab 


-15 K 


58 K 


-30 K 


-12 K 



TABLE IV. Normal mode coupling constants for the four 
Peierls-active phonon modes at qo, determined by Eq. (p2[). 



The solution of the equations gives the coupling constants 
as a function of N/(S cl0 ). The latter is then determined 



by the critical temperature Tsp 

(^qo) 



N 



N ^ 
l 



(-1) 



14.1 K via Eq. (|25|). 
Si • Si +f ) = 0.59 (33) 



For a spin- 1/2 system with two Cu chains per unit cell 
we have (Sq )/N < 0.75 where 0.75 is reached in the 
fully dimerized state. In the uniform Heisenberg case 

(Sq )/N = 0. 



In Table IV we show the calculated coupling constants 
of the spin system to the Peierls-active eigenmodes of 
the lattice at the wave vector of the instability qo- The 
signs are such that all contributions in the spin-phonon 
coupling term in the Hamiltonian ( p4| ) are negative when 
the phonon modes are macrosco pically occupied as de- 
termined via Eq. (H) in Sec. |VIIA| . The mode at 
£l2/(2n) = 6.5 THz is dominant, by its symmetry it es- 
sentially varies the angles r/. This will be reflected in the 
corresponding coupling constant discussed below. 

Note that the influence on the transition temperature 
of the lowest f2i/(27r) = 3.1 THz mode is as important 
as that of f2s/(27r) = 11 THz, because of the frequencies 
in the denominator of Eq. (|2£ 



VI. MICROSCOPIC COUPLING CONSTANTS 

The numerical values of the normal mode coupling con- 
stants thus given, the microscopic coupling constants can 
be determined. Using the matrix (^|) we rewrite expres- 
sion (|2C|) for q = q as 



(9L\ 

gL 
\g v oJ 



( 9iVth\ 

92 V^h 

93 V^h 



(34) 



and compute the coupling to the linear atomic elonga- 
tions. Then we calculate the angular and bond length 
couplings using Eq. fl23|). The resulting values are given 
in Tables and VI, respectively. 

The results allow for some immediate conclusions: 

• The coupling to 77, i.e., g r) , is the dominant contribu- 
tion. 

• The signs of the coupling constants are correct, J 
increases with increasing angles and decreasing 02-Ge 
bond length. The positive value of g£ a indicates that 
the ferromagnetic exchange is weakened more than the 



G 





aL/ k B 


go2/k B 


floa/feB 


Method [Reference] 


dJ 


dJ 


-890 K/A 


-110 K/A 


400 K/A 


-91 K/A 


Harm, theory and mean-field [here] 


10 %dk 





TABLE V. Coupling constants for the linear atomic dis- 
placements calculated via Eq. (^) using the values for g\ 
from Table 0. 



gs/ks 



gL/kB 



15 K/deg 



1.5 K/deg 



180 K/A 



-96 K/A 



TABLE VI. Coupling constants for the angles and bond 
lengths calculated via Eq. (^3|) using the values for g" from 
Table fv|. 



Microscopic superexchange fij] 22% 
Pressure vs. magnetostriction |l7| ~ 5% 



J dcg 



dcg 



dcg 



J deg 



- dcg 



TABLE VII. Variation of J with the variation of the an- 
gles. Note that in our notation there are two angles rf and 
two angles S K contributing each to the Cu-Cu superexchange 
path (see Fig. w). 



VII. STATIC DISTORTION 



antiferromagnetic exchange when stretching the 02-Cu 
bond. This is consistent with the net ferromagnetic ex- 
change of the 02-Cu plaquettes without the germanium 
side group predicted by Geertsma and Khomskii.li-3 

• Variation of the coupling constants shows g± to couple 
mainly to the angles <5, <?2 to u. and c/3 and 34 to be al- 
most entirely bond stretching.tj While the results for g v 
and gs are robust under variation of the parameters, the 
values of g^ u and g* e are less fixed within the accuracy 
of our approach. .— ■ 

• From magnetostriction data Buchner et al.il expect 
the influence of the Cu-02-Cu angle r\ on the magnetic 
exchange to be of the order of 2dJ/(Jdrj K ) « 10% 
per degree, and for the 02-02-Ge angle their value is 
2dJ/{Jd8 K ) « 1%. For J/k B = 150 K we obtain about 
twice the values (see Table |VIl| ) . 

• Comparing gf, B and <7go shows the effect of the ger- 
manium elongation on the magnetic exchange to be due 
mainly to the stretching of the 02-Ge bond. The con- 
tribution of the Ge side group to the magnetic exchange 
should depend on the 02-02-Ge angle as J sidc ~ cos 8. 
Therefore, the angle 8 K w 160° = 0.897T being close to ir, 
the angular dependency of J on 8 is quite small in spite 
of the large entire side-group effect, which is of similar, 
magnitude as that of the CUO4 plaquette elongation£M£l 

• Two groups analyzed the structural dependence of 
the superexchange within similar microscopic models. 
Geertsma and Khomskiili-3 obtained J/k B = 135 K and 
found 2dJ gco /{Jdif) « 16% and 2 <9J gco /( Jd8 K ) w 0.6% 
per degree. These values only accout for the "geometri- 
cal—, contribution and are thus lower bounds. Braden et 
alii found J/k B = 160 K and gave 2dJ/(Jdr] K ) « 44% 
and 2dJ/(JdS K ) 1.1% per degree. The agreement be- 
tween the microscopic models is affected by the choice 
of the parameters and the number of orbitals taken into 
consideration. 

A summary of the value s ob tained in the different ap- 
proaches is given in Table VII. 



The microscopic coupling constants given, we can di- 
rectly calculate the effect of static distortions of the lat- 
tice geometry on the magnetic exchange. 



A. Dimerization 

Using the static displacements of the ions in the spin- 
Peierls phase at T = 4 K,E-3 one may calculate the alter- 
nation of the magnetic exchange usually used in mean- 
field approaches to the spin-phonon coupling, i.e., 



J £ (l + (_!)'.*.,) S,,-S 



(35) 



This is achieved by substituting in the spin-phonon cou- 
pling term ( O ) the atomic displacements by their expec- 
tation values r°i — > (— l)'=+ /a = (t")t=4k and comparing 
the resulting (H™) T=4K with Eq. @. Equivalently 
one can calculate (Hf^} T= A^. by using the static angular 
and bond length deviational3 yielding the same results. 

We find 8 j J/k B = 17 K or Sj 0.11 (J/k B = 150 K). 
By solving the system of linear equations defined by Eq. 
( Pq ) for q = qo, the expectation values (b\)/VN have 
been determined to be 0.061, -0.11, 0.034, and 0.006 for 
A = 1, 2, 3, and 4, respectively. The elastic energy per 
unit cell of the spin-Peierls distortion at T ~ 4 K then is 
given by 



(Hp) 
Nk B 



E 



Nk T 



<6a) 2 = 5K. 



(36) 



This energy loss has to be compensated by the spin sys- 
tem. Considering that the maximum gain of magnetic en- 
ergy is reached in the fully dimerized case with 0.375 SjJ 
per Cu ion, we find a lower boundary for the dimeriza- 
tion of Sj > 0.044. Including a NNN term uf-jEq- @ 
with J2/J = 0.241 as studied by Chitra et alx3 using a 
density matrix renormalization group (DMRG) approach 
we find 8j > 0.078. 

Our result is within a factor of two of the values 
obtained by using dJ/dr] -apd dJ/dS obtained from 
the magnetostriction resultaiJ and from the microscopic 
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Method [Reference] 



5j 



Harmonic theory and mean-field [here] 0.11 

Macroscopic occupation of modes [here] >0.04 

Microscopic superexchange Jl4|,[l5| 0.07 to 0.2 

Dynamic phonons and experimental gap (l7| ~ 0.05 

Static phonons and experimental gap JTT| , p~2| 0.01 to 0.03 

Coupled chains juj 0.01 to 0.12 

TABLE VIII. Exchange-alternation in J[l + {-l) l '8j]. 

models. 00 All other published estimates of the dimer- 
ization result from an analysis of the magnetic excita- 
tion spectra observed by inelastic neutron or Raman 
scattering. Most of these estimates are based on the 
static dimerized Hamiltonian ([35]) (Refs. [ll]|lj) and yield 
dimerization va lues m uch smaller than the aae reported 
here (see Table, VIII ) . Augier and Poilblanccj as well as 
Wellein et alEU extend the static model by coupling to 
dynamical phonons which reduces the magnetic gap by 
lowerins-the effective lattice distortion acting on the spin 
system.tJ The derivation of their model and the signifi- 
cance of the phonon dynamics are more closely discussed 
in Sec. VIII. Introducing interchain coupling may further 
suppress the spin gap.Ej For an extensive discussion, see 
Ref. 

All methods incorporate more or less crude approxima- 
tions to the real physical situation leaving the question 
of the true value of Sj unanswered. Our lower boundary 
should be rather reliable though. The values obtained 



in the different approaches are given in Table VIII for 
comparison. 



B. Pressure 

Brauninger et aZ.S and Braden et aZ.0 have investi- 
gated the pressure dependence of the angles and bond 
lengths in CuGeC>3 under hydrostatic pressure. The lin- 
earity of the pressure dependence is reasonable for pres- 
sures < 2 GPa. The values for the pressure gradients 



obtained from Ref. 20 are shown in Table [X 



Regarding the partial derivative of the exchange in- 
tegral dJg/dp = (dJ/d9)(dO/dp) we find immediately 
the pressure gradients of the different angular and bond 



length contributions to J as given in Table IX 



Considering all four contributions, we obtain the total 
variation of the antiferromagnetic exchange. 



d J „<9J„ „dJ$ „<9J Go ,<9J Cu 
— = 2 — 2 + 2 — - + 2 — — + 4 — — 
dp dp dp dp dp 



= -9 



fc B K 
GPa 



(37) 



For J/k-Q = 150 K this value corresponds to 
dJ/(Jdp) w —6% per GPa. The pressure dependency 
of the magnetic susceptibility is directly related to the 
magnetostriction. A value of —dx/{x®P) ~ dJ/(Jdp) w 



dp 

9J e 
dp 



-0-16 ^ 



-2.5 



1 o dcg 
~ L -° GPa 



-1.9 



-0.0012 
-0.22 



GPa 



-0.0033 
0.32 



TABLE IX. Linear pressure gradients of angles and bond 
lengths from experimental data in Ref. ^ (top) and the result- 
ing theoretical pressure gradients dJg/dp — (dJ/d8)(d8/dp) 
(bottom). The values for dJ/89 are given in Table N& 



Method [Reference] 



3J 
Jdp 



Harmonic theory and mean-field [here] 
Susceptibility via magnetostriction || 
Curie- Weiss fit to the susceptibility Q 
Fit to the triplet dispersion flOl 



-6% 
-8% 
-7% 



GPa 
1 

GPa 
1 

GPa 



-8 to -10% ^ 



TABLE X. Variation of J with pressure. 



—8% per GPa was_r>btained after averaging the uni- 
axial components.aLj Takahashi et al.u have measured 
the pressure dependence of the Curie constant C by fit- 
ting a Curie- Weiss law to the high-temperature tail of the 
magnetic susceptibility. Assuming C ~ 1/ J one can esti- 
mate a value oLabout dJ/(Jdp) as —7% per GPa. Nishi 
and co-workeralj compared fits to the dispersion of the 
lowest triplet excitations at different pressures. They as- 
sume the ratio between the exchange J and next-nearest- 
neighbor exchange J 2 with a value of J2/J ~ 0.25 not to 
alter under pressure, and find dJ/(Jdp)r-y —10% per 
GPa. In contrast to that Fabricius et alEQ find J2 not 
to alter under pressure. Then the result from Nishi's 
analysis is corrected to dJ/ (Jdp) « —8% per GPa. A 
summary of the values is given in Table [x| showing their 
consistency. 



C. Thermal expansion and spontaneous strain 

In a harmonic lattice the coefficients of linear thermal 
expansion a = (dL) / (LdT) p vanish. Here L is the length 
of the crystal in a given spatial direction. Anharmonic 
contributions result in temperature-dependent phonon 
frequencies which in turn yield finite values for a. The 
coefficient of thermal expansion is linked to the specific 
heat via the (temperature dependent) Griineisen param- 
eter. This implies in the limiting cases T — > : a ~ T 3 
and T ^> Bb-,: a ~ constant, where Grj is the Debye 
temperatureH 

The thermal expansion in CuGe03 can be attributed 
to two effects: the usual anharmonic behavior de- 
scribed above and anomalies due to the spin-phonon 
couplingHj'Ea The coefficient of linear thermal expansion 
of the c axis in CuGe03 has a negative sign between Tsp 
and T ~ 200 K. The expansion of the c axis enlarges 
J via the angle 77. The spin system then gains energy 
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e 


V 




5 




ae K 

8T 


n 9 dc s 

250K 


0.6 


dog 
2. r ,()K 


0.0002 ^ 


8J e 
dT 


°' ± 250K 


0.9 


fe B K 
2 r > ) ! \ 


0-04^ 



dc 

-0.002 

0.2 



Cll C-2 



250K 



TABLE XI. Experimental linear temperature gradients of 
angles and bond lengths from Ref. ^ (top) and the resulting 
theoretical contributions to the temperature dependence of J 
(bottom) between 20 and 295 K. 



when the temperature is lowered to T ~ J by driving 
the anomaly.Efl A rough quantitative estimation can be 
extracted from the analysis of the temperature depen- 
dence of the herein copsidered bond lengths and angles 
given by Braden et alLn Their temperature dependence 
between 295 and 20 K is close to linear and presented in 
Table K3L Summing up the different contributions equiv- 
alently to Eq. @ yields d,J/(JdT) w -2.6% per 250 K 
(J/ks = 150 K). This effect is a superposition of the nor- 
mal thermal expansion with positive <9J norm / (JdT) and 
the anomalous effect at low temperature which can be 
estimated by dJ an /(JdT) < 2dJ r) /(JdT) = -4.1% per 
250 K. 

As the crystal undergoes the SP transition spon- 
taneous s train appears along all three orthorhombic 
directions OEJ The strain couples different T^ modest^ 
and gives a correction to Eq. (|32|) which we now show to 
be unimportant. The elastic energy per unit cell related 
to the spontaneous strain at T ~ 4 K can be estimated 
from the elastic constants. The diagonal elastic constants 
were taken from the ultrasound study by Poirier et aZ.,E3 
and off-diagonal rte rms were calculat ed w ith the lattice- 
dynamicaLmodelcJ as shown in Table XII using standard 
notation]^! With .the values for the strain e.j given by 
Winkelmann et alc3 we find 



k B 



i ■ b ■ • 
~2k^ 



E 

i,j=l,2,3 



etCijCj = 7 - lO- 



rC. 



(38) 



Note that the strain components e± = £5 = eg vanish, 
since the orthorhombicity is conserved. The elastic en- 
ergy involved in the strain is four orders of magnitude 
smaller than the elastic energy of the dimerization given 
in Eq. ©. 

Note that the components of the spontaneous straino 
have the opposite sign.-compared with the anomalies of 
the thermal expansion^] discussed above. The sponta- 
neous strain may thus be interpreted as a relaxation of 
the latter when the spin system changes its character at 
the spin-Peierls transition. The relaxation is of the order 
of 1%. 



VIII. COUPLING CONSTANTS FOR 
REAL-SPACE NORMAL COORDINATES 

In order to obtain real-space expressions we use the 
Fourier representation of the Bose operators 



Expt. [lO 11 ^ 

1 L cm 

Theory [lO 11 



dyn " 



7.4 
8.2 



2.1 
5.0 



C33 
33.2 
34.6 



Ci 



C, 



C 2 



3.0 4.0 2.2 



TABLE XII. Experimental uniaxial elastic constants from 
Ref. |37] (top) and the theoretical elastic constants obtained 
from the shell model (bottom). 



A 



4 



g^/ka -0.5 K 17 K -20 K 1.6 K 

g l r/kn -7.2 K 12 K 4.8 K -7.5 K 

TABLE XIII. Coupling constants to real-space normal 
modes obtained from Eqs. ( fio] ) and (p4|). 



(39) 



For simplicity we neglect the wave- vector dependence of 
the frequencies Q\ = Qx.q,, and of the polarization vec- 
tors e"fA) = e"(A, qo). The coupling constants g\(q_o) in 
Eq. (EG) then are divided into 



2VL> 



-9\ 



eS u (A) 



(40) 



Transforming the Hamiltonian (|24|) via Eq. (|3^) we ob- 
tain in real space 



A.l 



H r = J2 ( b{ tl b XA + - 

•Ec- 1 ) 1 ' 

A,l 

9? (b[, 



J ^2 Si-Si+i 



9? W 



,1 -r 6 A,1 



+ h \,\+y + h \,i+y 



Si • S 



(42) 



The coupling constants are given in Table Kill. This 
result implies that the oxygen and germanium displace- 
ments are of the same importance for the spin-phonon 
coupling as the copper elongation. 

Motivated by the symmetry of the Peierls-active 
phonon modes an effective one-dimensional model can 
be obtained by restricting the sum to a single chain. The 
Fourier transform of the one-dimensional model derived 
from Eq. ( |42"|) shows the different q dependences (q = q z ) 
of the copper and the local term. 



ff 1D = £ftn A [b{ q b x ,- q 



(43) 
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Here we defined the ID coupling constant 



9xM = {l-e^)9T + l9T- (44) 

Several studies@0iH have been carried out using 
real-space Hamiltonians in the form of Eq. ( f42| ) reduced 
to a one-dimensional model. Usually a single-mode 
Hamiltonian, only keeping the local term is considered, 
i.e., in their notation 2g l f c = g, while the other coupling 
constants are set to zero. Considering g% n , g 1 ^", and g^" 
being of the same order of magnitude, this simplification 
should only yield qualitative results. 

Yet, these treatments include the dynamics of the 
phonons. The significance of the latter can be estimated 
from the size of the zero-point motion of the ions. With- 
out the negligi ble con tribution from the macroscopic oc- 
cupation (Sec. VII A) the fluctuations of the modes 
at T = can be obtained from Eq. da) using the approx- 
imation of dispersionless phonons introduced above: 



((' 



E 

A 



e?(A) 



h 
2TV 



(45) 



The resulting values are J ((r%) 2 ) = 0.029, 0.035, 0.048, 

and 0.053 A for r" = r^, u , r% m , r* 2 , and r^, 2 , respectively. 
They are consistent with the values of the total zero- 
point fluctuations obtained from the shell model and the 
neutron-scattering experiments presented in Ref. |3^. The 
zero-point fluctuations are thus a factor of 5 to 10 larger 
than the static distortions as given in Eq. ([ll"|). 

On t he o ther hand, the Ginzburg criterion discussed 
in Sec. VA and the consistency of our results with ex- 
perimental ones justify our mean-field approach. In ac- 
cordance with that, Kliimper et alEB show that a vari- 
ety of physical quantities can be obtained correctly in a 
mean-field picture. It is beyond the scope of this paper 
but certainly an interesting question addressed to future 
studies which quantities are sensitive to the zero-point 
fluctuations and why. 



IX. DISCUSSION OF X o 



binding energy of the singlet baund state for T < Tsp, as 
seen by Raman experiments.Ej Interchain coupling will 
reduce the value of xo because.-af an enhancement of the 
antiferromagnetic correlations.E3 

As can be seen from Eq. ( p5[ ) our coupling constants 
scale as g\ ~ ^/Xo -1 - From the above results follows 

1 < -\ZXo~ 1 < an d we adapt the mean value of 
Xo = 0.5 for our calculations. This value is close to 
the result for J^/J = 0.241. Within the accuracy of 
our approach we can use J = 150 K as given by Castilla 
et al.iB The choice of xo is justified a posteriori by the 
agreement of the results in the literature. Also note that 
including a NNN term with J 2 /J = 0.24 in Eq. @ with 
5 j = 0.1 and using exact diagonalization gives a value 
of E; a ( s h ■ s h+i)/ N = 0.57 per two Cu ions in 

agreement with the value given in Eq. (|33]). 

Applying hydrostatic pressure the transition temper- 
ature grows at a rate of 4.8 K/GPa.B In our approach 
Tsp is given by Eq. ( p5| ) and depends on the coupling 
constants g\, the frequencies fix, and the factor xo- The 
coupling constants g\ in turn depend on the linear deriva- 
tives of the magnetic exchange and the polarization 
vectors, as given in Eq. (|2(]). In a harmonic lattice the 
phonon frequencies and polarization vectors are indepen- 
dent of pressure. It seems very unlikely that the Peierls- 
active modes exhibit extremely large negative Griineisen 
Parameters which would be needed in order to describe 
the increase of Tsp upon pressure via the pressure de- 
pendence of the phonon frequencies. The linear coupling 
constants g™ also are independent of pressure, and since 
the lattice distortions are rather small,E3E3 we do not 
expect higher-order contributions to play a crucial role. 
We must thus conclude the value of xo to be strongly 
pressure dependent. 

Together with the pressure dependence of J2/J dis- 
cussed by Fabricius et aL,E3 this may explain the shift of 
Tsp .Ej When introducing interchain coupling, prefactors 
and the functional dependence-pf the spin-spin correla- 
tion function are also alteredo The compressibility of 
the crystal is largest in the b direction so that the alterna- 
tion of the interchain coupling under pressure is another 
possible origin of the pressure dependence of xo and Tsp. 



The approach by Cross and Fisheruu gave a value of 
Xo ~ 0.26. This value is independent of J because of 
the scale invariance at q c = n/c. The scaling hypoth- 
esis is applicable close to the critical point of the spin 
chain, i.e., in the limit T^a-Q. Recent DMRG results ob- 
tained by Kliimper et alr3t?3 show a strong temperature 
dependence of Xo(T/ J). For J 2 = and J = 120 K they 
found Xo(T S p/J) w 0.28. For J 2 /J = 0.241 and J = 150 
K the parameter attains Xo(Tsp/J) ~ 0.56, whereas for 
J 2 /J = 0.35 and J = 160 K they found Xo(T S p/J) « 1. 

The exact value of J2 in CuGe03 has not yet been de- 
termined. Fits tOpthe-susceptibility for T > Tsp indicate 
an overcritical J2 jii'EII but fits tenthe four-spinon contin- 
uum seen by Raman scattering^ indicate an undercrit- 
ical J2. In favor of an undercritical J 2 is also the small 



X. SUMMARY 

In this paper we have given a detailed analysis of the 
microscopic magnetoelastic coupling in CuGeOa which 
may be easily transferred to other systems. The compar- 
ison of several theoretical and experimental approaches 
yields a satisfa ctory consistency. Numbers have been 
given in Table VII for the a ngula r dependence of the 
magnetic exchange, in Table VIII for the dimerization, 
and in Table [x] for the pressure dependence of the mag- 
netic exchange. The quantitative agreement of course is 
limited by the uncertainties within experiments and the- 
oretical techniques. Coupling constants for effective one- 
dimensional real-space model Hamiltonians accessible to 
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numerical studies are given in Table^H W e have dis- 
cusse d the applicability of static models (Sec. VII A and 
VIII), and we were able to explain qualit atively the c-axis 
anomaly of the thermal expansion (Sec. VII C ) . 
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